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Разработана физико-математическая модель массопереноса при ионно-лучевой кристаллизации. Рассмотре- 
но получение фотоактивных структур $(п’ /5Кр)/5Кр*) на подложках 100 мм методом ионно-лучевой кри- 
сталлизации. Определены оптимальные условия этого процесса: остаточное давление в ростовой камере — 
10`* Па; температура подложки — 550 °С; плотность ионного тока — 2 мА/см’; ускоряющее напряжение пуч- 
ка — 400 В; расстояние «мишень — подложка» — 150 мм. Данные сканирующей электронной микроскопии 
показывают, что выращенные фотоструктуры имеют практически бездефектную поверхность. Установлено, 
что фотоструктуры 5(п* У5КрУ$Кр”) характеризуются внешним квантовым выходом более 90 % в диапазоне 
длин волн 550—900 нм (спектр АМ 1,5) при следующих условиях: толщина фронтального слоя 5(п*) 100 нм, 
уровень легирования п* = 5'10' см 3, толщина слоя 5(р) 130 мкм, уровень легирования р = 2'10'° см-3, 
толщина слоя 5(р*) 500 нм, уровень легирования р* = 1'108 см 3. 

Ключевые слова: ионно-лучевая кристаллизация, фотоактивная структура, внешний квантовый выход. 


Введение. Достижения в фотовольтаике стимулировали исследования фотоактивных материалов 
и структур, а также разработку методов их получения. На основе кремния и соединений АТВ" со- 
зданы высокоэффективные солнечные элементы [1]. 

Для получения структур с фотоактивными областями широко применяются методы моле- 
кулярно-лучевой [2] и газофазной эпитаксии [3]. В последние три десятилетия технологии полу- 
чения фотоактивных материалов и структур стремительно развиваются, производство становится 
сложнее. При этом исследователи продолжают искать новые методы получения фотоактивных 
материалов. Цель настоящей работы — получение и исследование фотоактивных слоёв и струк- 
тур на основе $ методом ионно-лучевой кристаллизации (ИЛК). 

Теоретическая часть. Для теоретического исследования ионно-лучевой кристаллизации одно- 
компонентных полупроводниковых материалов разработана математическая модель процесса. 
Она основана на имитационном подходе Монте-Карло. 

На рис. 1 показана схема моделируемого процесса ионно-лучевой кристаллизации и ис- 

пользуемые математические величины. Первичный пучок ионов аргона с энергией Е’. , плотно- 


стью тока ],., диаметром @. падает на центральную часть мишени. Причём диаметр мишени 


Аг* ! 


должен быть больше диаметра ионного пучка. Угол падения (между вектором /. и плоскостью 


мишени) — произвольный а = (20—70)°. Пучок ионов выбивает атомы мишени, центр которой 


* Работа выполнена при поддержке Министерства образования и науки РФ в рамках ФЦП «Исследования и разработки по 
приоритетным направлениям развития научно-технологического комплекса России» на 2007—2013 годы (ГК 
№ 14.516.11.0012) и ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009—2013 годы (грант 
№ 14.В37.21.1105). 
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удалён от подложки диаметром Д.„ь (с радиусом К») на расстояние [. Количество атомов росто- 


вого материала определяется коэффициентом распыления, зависящим от Е, иа. 
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Рис. 1. Схема процесса ионно-лучевой кристаллизации 


Условия и симметрия задачи допускают одномерное моделирование процесса при следу- 
ющих допущениях: 
1) распылённые атомы распределяются в пространстве по закону косинуса; 
2) распылённые атомы не сталкиваются друг с другом и с атомами рабочего газа; 
3) распылённые атомы конденсируются в точке соударения с подложкой; 
4) максимальная интенсивность выбитых атомов наблюдается при углах распыления, близких 
(или равных) углу зеркального отражения; 
5) ионы первичного пучка имеют одинаковую энергию и распределяются по сечению ионного 
пучка равномерно. 


Количество ионов /М„ с зарядом 9 в ионном пучке сечением 5 и плотностью тока },., 


провзаимодействовавших с поверхностью мишени за время проведения процесса т, устанавли- 
вается выражением: 


от 
лы (1) 


Предполагается, что ионы в пучке движутся по равноудалённым параллельным траекториям и 
сосредоточены в р сечениях. Число ионов в одном сечении равно: 
5 — М: (2) 
р 
Процесс моделирования является итерационным. Так как ионы первичного пучка моноэнергетич- 


ны и распределены по его сечению равномерно, то для получения координат локальной области 
распыления атомов мишени можно воспользоваться генератором псевдослучайных чисел, подчи- 
няющимся равновероятностному закону распределения. Координата точки распыления Ёго атома 


Гоп 


х„„/ выбирается на отрезке [#„., Ё „м ] при помощи генератора равновероятностного распреде- 
ления: 

Хриь г = (Е = | -“,, (3) 
где /=1,№5 , , — случайное число, равномерно распределённое в интервале [0;1]. 


Далее для каждого распылённого атома рассчитывается траектория движения в системе 
отсчёта ху (рис. 1, а). Угол, задающий направление движения в системе ху, моделируется генера- 
тором псевдослучайных чисел, подчиняющимся косинусоидальному закону распределения: 

9.„„, = агссо$ (2%, 1), (4) 


5$риё! 
где /=1,№°. 


оп 
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Уравнение траектории атома, распылённого из точки с координатами (х 0), имеет вид: 


5$риЁ! ! 


я. =19(6.„,,)(х + Хи). (5) 
Уравнение подложки в системе ху имеет вид: 
ГВ 
со$ (а) ` 
Возможная координата кристаллизации атома на подложке х 


у =®9(а)х + (6) 


находится решением системы 


сопа ! 


уравнений, образованной выражениями (5) и (6): 


ЕЙ (7) 
сопа ! к, те К, 
где К, = 9(а), = {9 (агссоз (2%, -1)), Ь Ра и, — случайное число, равномерно распре- 


делённое в интервале [0; 1]. 
Если х Е ра Жо в] ‚ то атом считается осаждённым на подложке. В противном 


сопа ! 


случае атом не взаимодействует с подложкой и не учитывается при моделировании роста слоя. 

Для построения профилей распределения толщины выращиваемых слоёв подложка раз- 
бивалась на Л одинаковых ячеек длиной /м. Номер ячейки / в которую конденсировался атом, 
определялся из выражения: 


а ВЕ ЖЕ 


_х 


=Е| п 





, (8) 
(х 506 _ ПОВЁ 5иь _ ей ) 
где Е — функция округления. 

При повторении описанной процедуры для № „ подсчитывалось общее число распылённых 


атомов мишени //.„,, конденсированных атомов №, число атомов в каждой ячейке №», 


определялись коэффициент массопереноса (9) и толщина слоя в /й ячейке (10): 
М 


Ки =, (9) 
Е 
м, КкК.у..т 
р, = сопа! т а а т | (10) 
и р. 9 


где т. и ра — атомная масса и плотность распыляемого материала. 


с-51(и”) — 0,1 мкм 


с-5Ир) 130 мкм 


с-З (р = 0,5 мкм 


Рис. 2. Выращиваемая фотоактивная структура 





На рис. 2 представлено схематическое изображение выращиваемых фотоактивных струк- 
тур. В работе [4] проведено компьютерное моделирование фотоактивных структур предложенной 
конструкции и показано, что эффективность выращиваемых фотоактивных структур будет макси- 
мальной при толщине фронтального фотоактивного слоя 9(л”) 100-300 нм и уровне концентра- 
ции легирующей примеси 5'10*8 см-3. 
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Экспериментальная часть. Для исследования свойств фотоактивных структур 5(л”)/5(р)/5(р") 
методом ионно-лучевой кристаллизации были изготовлены три серии образцов. Процессы прохо- 
дили при различных температурах подложки, энергиях распыляющих ионов и временах после 
диффузионного отжига. В качестве основы для базовой области структуры $(л“)/5(р)/5(р") были 
выбраны подложки кремния р-типа проводимости с ориентацией (100). Мишенями для выращива- 
ния фронтального слоя служили пластины низкоомного кремния /л-типа проводимости с ориента- 
цией (100). 

Формирование высоколегированных областей л’ и р* проводилось стандартным методом 
диффузии из легированного слоя $0.. Получение р"-слоя структуры осуществлялось при темпе- 
ратуре 950 °С в течение 40 минут в атмосфере кислорода. Фронтальный фотоактивный слой 
т-типа формировался при температуре 950 °С в течение 10 минут в атмосфере кислорода. Глу- 
бина залегания р-п-переходов контролировалась измерением вольт-фарадных характеристик ба- 
рьера Шоттки, образованного ртутным зондом и полупроводниковой структурой. Структура и 
морфология поверхности исследовались методами сканирующей зондовой микроскопии (СЗМ). 
Фотоэлектрические свойства структур исследовались измерением спектральных зависимостей 
внешнего квантового выхода. 

Результаты экспериментов и их обсуждение. Ранее в работе [5] были исследованы зависи- 
мости процесса ИЛК от параметров конфигурации системы «источник ионов — мишень — под- 
ложка». Установлено, что слои с минимальным радиальным разбросом толщины будут осаждаться 
на подложки 100 мм при следующих условиях: 

— плотность ионного тока — 2 мА/см?; 

— расстояние «мишень — подложка» — 150 мм; 

— угол наклона «мишень — ионный пучок» — 45°. 

Методами сканирующей зондовой микроскопии исследовано влияние температуры под- 
ложки и энергии ионов пучка на структуру и качество выращенных методом ИЛК слоёв кремния. 
Видно, что при температуре 550 °С (энергия ионов фиксирована и равна Е „„‚ = 380 эВ) получа- 


ются практически бездефектные и однородные слои кремния высокого качества (рис. 3). Это 
можно объяснить следующим образом. 

При ионно-лучевой кристаллизации для создания условий равномерного роста слоя необ- 
ходимо обеспечить достаточную скорость осаждения материала, чтобы при данной температуре 
происходил равномерный послойный рост плёнки на всей площади подложки $1. На основе дан- 
ных (рис. 3, а, 6, в) можно сделать вывод, что при температуре подложки 550 °С (рис. 3, 6) реа- 
лизуется оптимальное соотношение скорости поверхностной диффузии атомов к скорости оса- 
ждения распылённых атомов. 

Результаты исследования влияния энергии ионов на качество формируемых кремниевых 
слоёв при температуре подложки 550 °С отражают СЗМ-изображения поверхности центральной 
части слоёв при трёх различных значениях энергии распыляющих ионов (рис. 4, а, 6, в). При 
энергии ионов 400 эВ наблюдается практически эпитаксиальный рост малодефектных слоёв 
кремния (рис. 4, 6). При энергиях ионов, отличных от 400 эВ, превалируют процессы дефектооб- 
разования и ухудшения структуры выращенного слоя (рис. 4, а, в). 

Проведены СЗМ-исследования влияния температуры подложки и энергии ионов на струк- 
туру поверхности осаждённых слоёв кремния. Данные этих экспериментов позволили установить 
оптимальные параметры ростового режима и конфигурации системы «источник ионов — ми- 
шень — подложка», при которых выращенные слои кремния обладают высоким структурным со- 
вершенством и минимальным разбросом толщины: 

— давление в ростовой камере — 10-4 Па; 
— температура подложки — 550 °С; 
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— плотность ионного тока — 2 мА/см’; 

— энергия ионов пучка — 400 эВ; 

— угол наклона «мишень — ионный пучок» — 45°; 
— расстояние «мишень — подложка» — 150 мм; 
— диаметр подложки — 100 мм; 

— диаметр ионного пучка — 44 мм. 


= 





‚и 6.04 _._ 
З мкм 


%® = Е ЕЕЕЕЕЕЕЕЕЕНЕ 


5 мкм 
в) в) 
Рис. 3. Результаты СЗМ-исследований влияния температуры Рис. 4. Результаты СЗМ-исследований влияния энергии ионов 
на качество осаждаемых слоёв (энергия ионов Е’. = 380 ЭВ ): на качество осаждаемых слоёв (температура 7 = 550 °С): 


Т= 400°С (а), Т= 550 °С (6), Т= 700 °С (8 Е‚, =150 38 (а), Е’. =4003В (6), Е. =600 3В (в) 
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Рис. 5. Результаты исследования структуры 5(л*)/5(р)/5Кр\): глубина залегания р-п-перехода (а), спектральная зависи- 
мость внешнего квантового выхода (6) 


Исследован процесс формирования р-п-переходов методом ионно-лучевой кристаллиза- 
ции для получения фотоактивных структур. Результаты представлены на рис. 5, а. Анализ пока- 
зывает, что увеличение времени отжига приводит к улучшению равномерности распределения 
фосфора по глубине фронтального слоя 5(л*). Однако с улучшением однородности распределе- 
ния концентрации фосфора по глубине фронтального слоя наблюдается смещение границы р-и- 
перехода в объём структуры, сопровождающееся уменьшением концентрации фосфора в фрон- 
тальном слое. Экспериментально установлено, что при отжиге в течение 10 минут наблюдается 
оптимальная глубина залегания р-п-перехода (100—200 нм) в сочетании с высоким уровнем леги- 
рования (4`10'8 см-3) и равномерным распределением концентрации фосфора по глубине фрон- 
тального слоя 5’). 

Для оценки качества и анализа эффективности фотоэлектрического преобразования излу- 
чения в фотоактивных структурах 5(')/5(р)/5К( р"), полученных ионно-лучевой кристаллизацией, 
выполнены исследования спектральных зависимостей внешнего квантового выхода (рис. 5, 6). 
Полученная спектральная зависимость характеризуется высоким внешним квантовым выходом в 
спектральном диапазоне 550—900 нм, который достигает значений до 90 %. Это, по-видимому, 
обусловлено высоким структурным совершенством фронтальных слоёв $, выращенных ионно- 
лучевой кристаллизацией. Полученные экспериментальные данные хорошо согласуются с резуль- 
татами компьютерного моделирования [4] фотоактивных структур 5$(*)/5(р)/5Кр”) для толщины 
фронтального слоя 5(л”) 0,1 мкм и уровня легирования 51018 см-3 в спектральном диапазоне 
300—1200 нм при условии освещения АМ 1,5 (рис. 5, 6 — сплошная линия). 

Выводы. Разработана математическая модель ионно-лучевой кристаллизации фотоактивных ма- 
териалов, которая позволила теоретически исследовать зависимость процесса ионно-лучевой 
кристаллизации слоёв от параметров конфигурации системы «источник ионов — мишень — под- 
ложка». На основе теоретических данных экспериментально определены оптимальные параметры 
формирования слоёв $ для получения фотоактивных структур 5(л“”)/5(р)/5(р”) методом ионно- 
лучевой кристаллизации. Установлено, что однородные по толщине, практически бездефектные 
слои кремния на подложках 100 мм получаются при следующих ростовых режимах: 

— давление в ростовой камере — 10-4 Па; 

— температура подложки — 550 °С; 

— плотность ионного тока — 2 мА/см?; 

— ускоряющее напряжение пучка — 400 В; 
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— угол наклона «мишень — ионный пучок» — 45°; 
— расстояние «мишень — подложка» — 150 мм. 

Данные анализа профилей распределения концентрации фосфора в фронтальном слое с- 
9(л’) свидетельствуют о том, что увеличение времени отжига приводит к улучшению равномер- 
ности распределения фосфора по глубине фронтального слоя с-$(л”). Экспериментально уста- 
новлено, что оптимальная глубина залегания р-п-перехода (0,1—0,2 мкм) наблюдается при отжи- 
ге в течение 10 минут при температуре 900 °С. 
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ОЕВТУАТТ2АТТОМ ОЕ РНОТОАСТТ\Е $ТВУСТУВЕС $1( №) /$Т(Р)/$Т(Р!) ТНВОУСН 
ТОМ-ВЕАМ СВУЗТАИТ?АТТОМ" 


\У. М. Уагаука 

(Боп За Теспптса! Упм№ег$Ку), 

$. М. СпеБокагем, А. $. РазпспепкКо 

(бощиегпт баепИЯс Сепге ог ВизЗап Асадету оЁ 5аепсеэ), 
\У. А. ТКПа 

(бреда! Вечап-Епдтеенпа Вигеаи “Мпуегяуа” (С) 


Тре рпузса! ап та етайса! таз5-Напзег тодЕ! оЁ пе гоп-Беат сгуяаИРаНоп тео 15 аеуеореа. Тве 
депмайгаНоп о! рпоюасвуе ис шигез 5(п* )/5Кр/5р*) оп 100 тт зирягае$ Егоидй юп-Беат сгузеа!гайоп 15 
сопз/аегеа, Тве оритит сопа#опз Юг #е ргосезс: геядиа! ргеззиге п те дгоийй србатЬег — 10“ Ра; те 
зибзнае етрегашге — 550 °С; (юп ситепЁ аеп$йу — 2 тА/спт?; ассеегавоп уо№аде — 400 И: агдеЁ — 
5ир5иейе оапсе — 150 тт аге ащегттед. Тйе $сапптд еестоп тисгозсору ааа $Вои/ {паЁ {те дгоит 
рпоюасНуе ${тисигез Пауе созе а аиШе5$ 5итасе. ТПе гези!5 детопзв‘а аЁ те рпоюасвуе $гисЁигез 
$(п’/5КрУ$Кр”) ойег те ежета! диатит ейсепсу оуег 90 % т те иауеепдЁ гапде оЁ 550—900 пт 
(срейгит АМ 1.5) ипаег {те оюиита сопаоп5: Попе 1ауег И/сКпез$ 5(п*) — 100 пт; аопог аортда аепяйу 
П* = 5.108 ст 3; |ауег Искпез5 5(р) — 130 ит; ассерюг аорта аепэйу р = 2.10 ст; |вуег Нискпезя 
5Кр*) — 500 пт; ассерюг дорта аепзйу р* = 1.10: ст 3. 

Кеуигогаб: гоп-Беат сгузЕИ!гайоп, ррооасйуе Ягисге, ежегпа! диапит ейсепсу. 


* ТНе гезеагсН 15 допе ИН 1е зиррок тот е Митету оЁ ЕдисаНоп апа 54епсе оЁ пе Визчап ЕедегаНоп м/т {Не #ате оЁ 
{пе Редега! Тагдее Ргодгатз$ “В&0О т рйогКу деуеюртепЕ Яе!4$ ог {Пе э4епсе апа фесппооду эесог оЁ Кизча” Гог 2007—2013 
(Заке СопёгасЕ по. 14.516.11.0012) апа “Везеагсй ап асадетюс ха оЁ шпоуай\е Виза” Тог 2009—2013 (9дгапЕ 
по. 14.В37.21.1105). 


84 


